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Abstract 

We investigate the interface dynamics of the two-dimensional stochastic Ising 
model in an external field under helicoidal boundary conditions. At sufficiently 
low temperatures and fields, the dynamics of the interface is described by an 
exactly solvable high-spin asymmetric quantum Hamiltonian that is the in- 
finitesimal generator of the zero range process. Generally, the critical dynam- 
ics of the interface fluctuations is in the Kardar-Parisi-Zhang universality class 
of critical behavior. We remark that a whole family of RSOS interface mod- 
els similar to the Ising interface model investigated here can be described by 
exactly solvable restricted high-spin quantum XXZ-type Hamiltonians. 
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1 Introduction 



Mappings from two-dimensional (2D) Ising model interface configurations to diffusion pro- 
cesses are known at least since the work of Rost [1], and have been explored many times since 
then [2, 3, 4, 5, 6, 7, 8, 9, 10, 1 1]. In some cases [2, 3, 4], the interface dynamics of the model 
at zero temperature in the absence of an external field was mapped into the one-dimensional 
symmetric simple exclusion process, with the main result being the solution of a first-passage 
time problem showing that the mean lifetime of a shrinking domain is proportional to its initial 
area, providing a microscopic derivation for this well known experimental fact. It was also 
recognized that the resulting exclusion process can be recast as a probabilistic cellular automa- 
ton with a transition matrix equivalent to the transfer matrix (in a diagonal direction) of the 
symmetric six- vertex model in one of its critical lines [3, 4]. 

The relationship between interfaces, exclusion processes, and vertex models was explored 
further [12, 13, 14, 15], and it was realized that a Heisenberg Hamiltonian with pure imaginary 
Dzyaloshinsky-Moriya interaction that commutes with the transfer matrix of a six- vertex model 
describes the single-step surface growth model [16, 17], as well as a discrete-velocity version 
of the noisy Burgers equation, which in turn is equivalent to the Kardar-Parisi -Zhang equation 
[18]. In the interacting particle system scenario, the above mentioned Hamiltonian is but the 
infinitesimal generator of the asymmetric simple exclusion process [12, 14]. Conversely, a 
host of results concerning symmetric and asymmetric simple exclusion processes has been 
translated into the Ising interface problem and, in particular, the investigation of the motion of 
tagged particles, as first introduced in the study of the hydrodynamic behavior of exclusion- 
type processes [19, 20], has provided a partial explanation for the relationship between the 
characteristics of different asymptotic growth regimes in some (1 + 1) -dimensional stochastic 
growth models [6, 7, 9]. 

In this article we show that under suitable generalized, but otherwise quite natural peri- 
odic boundary conditions the dynamics of an interface in the 2D stochastic Ising model in 
the presence of an external driving field can be mapped via a particle-height transformation 
into the dynamics of hopping particles without exclusion known as the zero range process 
[21, 22, 23, 24, 25, 26, 27]. The infinitesimal generator of the zero range process is equiva- 
lent to a high-spin, in general asymmetric quantum Hamiltonian that is exactly solvable by the 
Bethe ansatz [28, 29, 30]. We argue that the critical behavior of a generalized particle-height 
model must be on the Kardar-Parisi-Zhang universality class of critical behavior, since this is 
the critical behavior of the corresponding generalized quantum chains. This may have implica- 
tions in the study of related models such as the dynamics of fc-mers and other Ising-type lattice 
configurations. 
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The article is organized as follows. In Section 2 we introduce the 2D Ising model in an 
external field and the single-spin flip rates in terms of which the dynamics of the Ising contours 
will be analyzed, and in Section 3 we show that it can be described by diffusing particles 
without exclusion and exhibit the infinitesimal generator of the process. Section 4 contains 
a brief exposition of the exact solution of the zero range process by the Bethe ansatz and a 
discussion on its dynamical critical exponent. In Section 5, we show that a whole class of 
interface models similar to the Ising interface model can be described by exactly solvable, 
generalized restricted XXZ-type Hamiltonians with many modeling possibilities. Finally, in 
Section 6 we summarize our results and indicate some directions for further investigation. 



2 The 2D stochastic Ising model in a field 

The 2D Ising model in an external field is described by the Hamiltonian 

H{S) = -J SrSr'-BY,Sr, (1) 

(r,r') r 

where S = {S'r : r G A^} with ^rGj— l,+l}are Ising spins, C is a member of a family 
of semi-infinite lattices of |A^| = L x oo sites, and (r,r') denotes pairs of nearest neighbor sites 
on A^. The integer index N in A^ refers to the boundary conditions, that are free in the 
infinite direction and helicoidal with pitch N in the finite direction, i.e., r + Lx + A^y = r for 
all r G A^. When N = we recover the usual periodic boundary condition, which is however 
uninteresting for our purposes, as we will see later. In the above Hamiltonian we take J > 0, 
making the model ferromagnetic, and for definiteness we take B ^ 0. 

We introduce a dynamics on the Ising spins through the master equation 



Sen(A~) 



W{S S)P{S,t)-W{S S)P(S,0 



(2) 



for the probability P{S,t) of observing the configuration S G ^^(A^) = { — 1,+!}'^^ at instant 
t, where W(S — S) is the rate at which configuration S is reached from configuration S per unit 
time. The rates VK(S — )• S) should be translation invariant and verify the condition of detailed 
balance W{S S)P(S) = W{S S)P(S), with P(S) oc exp[-j3//(S)] the Gibbs equilibrium 
probability distribution and where for the sake of notational economy we omitted the depen- 
dence of P{S) and W(S — )• S) on J, B, and the inverse temperature P = l/kBT.ln this work we 
consider heat-bath single-spin flip transition rates given by 

W(S ^ S) = W{Sr ^ Sr) = ^exp[-j8//(5r)], (3) 
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with 

H{Sr) = -J S,Sr'-BSr and Zr = ^ exp[-j8//(5r)], (4) 

(r':r) 5r=±l 

where (r' : r) = {r' € : |r' - r|= 1}. 

Let w{Sr) = 5 L(r':r) l-^r' " "^rl be the number of spins neighboring Sr that have the sign 
opposite to it. In terms of this quantity, the single-spin flip rates read 

W{Sr ^ Sr) = 1 +exp[4j87(2- w(5r)) + 2j8B5r] ' 

At sufficiently low temperatures, as long as B < 27 spins with w{Sr) = 0, 1 will hardly flip, 
because their transition rates become exponentially small when compared with the other rates, 
of the order of exp[— 2jS(27 — B)] at maximum. Processes with w{Sr) = 3,4 correspond to 
fast processes, since at sufficiently low temperatures and again in the range B <2J their rates 
become close to unity, W{Sr — )■ 5r) ^ 1 — exp[— 2j8(27 — B)]. In the low temperature limit and 
in the range B < 27, thus, the heat-bath single-spin flip rates (5) define a process in which only 
spins with w{Sr) ^ 2 have an appreciable flipping rate, and henceforth we ignore the flipping 
of spins with w{Sr) < 2.' Spins with w{Sr) = 3,4, in turn, can be avoided by choosing initial 
configurations in which the "+" phase is separated from the "— " phase by a single-valued, 
non-self-intersecting staircase-like interface as in Figure 1 . With initial configurations of this 
type and within the low temperatures and fields regime, we are left with a process in which only 
spins with w{Sr) = 2 flip. Since spins with w{Sr) = 2 lie at the interface, the above-defined spin 
flip dynamics actually defines an interface dynamics. In the next section we map this dynamics 
into an interacting particle system on the integers. 



3 Mapping to the zero range process 

Various possibilities exist to map the dynamics of 2D Ising interfaces into a system of inter- 
acting particles on the line. One possible map is obtained by associating with every vertical 
bond in the interface a particle and with every horizontal bond a hole [3]. In this way we end 
up with a simple exclusion process in which particles hop in an augmented one-dimensional 
lattice. Another possibility was given in [6, 7]. In this case, one considers a set of L particles 
on a one-dimensional lattice occupying the positions x^, I ^ £ ^ L. If one associates the par- 
ticle labels i with a horizontal coordinate, and the particle positions Xi with the heights of an 

^In realistic pseudo-two-dimensional S = ^ Ising-like materials, e.g. in the antiferromagnetic compounds 
K2XF4 with X = Mn, Fe, Co, or Ni, J/ks - 1-100 K, such that pJ > 1 implies T < IK [31]. The values for 
which B <2J thus lie in the range B < 1 .5 T, of the order of half the magnetic field strength of a typical medical 
MRI system. 
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Figure 1: Staircase-like 2D Ising interface separating the "+" and the "— " phases. In this figure, L = 18 
and = 15. When B > (B < 0), the "+" ("-") phase invades the ("+") phase, while for B = 
the interface only fluctuates about its initial shape. When j3 — > 0° and |B|< 27, the number of bonds 
is conserved. 

interface, then one has a one-to-one map between the set of particles and an interface. With the 
additional constraint x^+i —X(^l, the resulting model was called the particle-height model. 
The particle-height model thus establishes a map between the low-temperature dynamics of 
an Ising interface and the simple exclusion process, although the constraint on the particles 
positions seems a bit artificial in the interface scenario. 

Our mapping of the Ising interface to a set of hopping particles on the integers is equivalent 
to the particle-height mapping with Xf+i — ^ 0. Let be the set of all single- valued, non- 
self-intersecting staircase-like Ising interface configurations of length \ 'y^\ = L + N, N 0, 
in the infinite strip of width L with helicoidal boundary condition r + Lx + A^y = r. The 
dynamics of the interfaces under the action of the flipping rates W{Sr) in the regime where 
_). 00 and B < 2J preserves the length of the interfaces, i.e., given an initial configuration 
y^{t = 0) G r^, all subsequent configurations y^{t ^ 0) € F^. The number N of vertical 
bonds in the interface thus decomposes the state space of interface configurations in an infinite 
number of disjoint sectors, Fi = {JN^Qr^. 

The interface configurations are single-valued functions with respect to the horizontal 
coordinate. Let hi ^ I ^ i ^ L denote the height of the Ising interface at site £. Then the 
heights differences rii = hi^i —h^ are nonnegative due to the special form of y^, and their 
dynamics relates to the dynamics of the interface as follows. Let us consider the case < B < 
27, and that the "+" phase lies below the "— " phase, as in Figure 1. In this way, every "— " 
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("+") spin that flips contributes to the growth of the "+" ("— ") phase, increasing (decreasing) 
the height variable associated with its horizontal position by one unit. The single-spin flipping 
rates are given by 

When B = 0, we have a spin (time) reversal invariant system, and the interface does not move 
bodily. Otherwise, when < B <2J and j3 — )• oo, /? = 1 and ^ = 0, and the surface can only 
grow. The general situation 0<q<p<lis obtained by taking jS — )• oo and B — )• with 
j8B = constant. The heights differences moves corresponding to the flipping processes are 
{nf_i,nf) — — l,n£ + 1) when h( ^ h( + 1, and {n(_\,n() — )• + — 1) when 
hi ^ hi — I. The maximum possible height difference in is rii = N, in which case hi = 
h2 = . . ■ = hg = and hi^i = hi^2 = ■ ■■ = hi = N fox some I. More generally, we have 
Y!eZ\ ng = ht — hi = N. Now one appreciates the role of the helicoidal boundary conditions 
on A^: the pitch N gives the total number of interacting particles in the heights differences 
scenario. The boundary condition in the heights differences scenario is simply periodic. 

The dynamics of the variables n£ is but the zero range dynamics, in which particles hop on 
the lattice without exclusion [21, 22, 23, 24, 25, 26, 27]. Except for the (immaterial) absolute 
values of the heights, the dynamics of the rii variables contains all the information about the 
evolving Ising interface. The elementary processes for the nc variables are 

P 

{ni^i + \,ni) ^ {ne^i,ni + \), ^ ni^i,ni ^ N - \, 2 ^ i ^ L. (7) 
<? 

As is well known [33, 34], we may write the master equation for reaction-diffusion processes 
on the lattice as a Schrodinger-like equation in Euclidean time, the infinitesimal generator of 
the Markov semigroup playing the role of the quantum Hamiltonian. In this scenario, the 
infinitesimal generator of the above zero range process is given by the x A^^ matrix operator 

L N-l 



e=im+n=0 



(8) 



where S^'" = 1 (gi • • • (gi 1 (gi S™'" (g 1 (g • • • (g 1, with 1 the x identity matrix and {E'"''')ij = 
Si.m5j^„ the N X N matrix with a single unit element in row m and column n occupying the 
^-th position in the direct product. The operator Hf^ can be diagonalized by the coordinate 
Bethe ansatz and, indeed, it has been diagonalized by this and related methods many times in 
the literature [28, 29, 30, 32]. We will thus not reproduce a complete resolution of (8) here. 
Instead, we just outline the technique and quote the main results regarding the process defined 
by Hn of interest to us. 
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4 Bethe ansatz solution 



4.1 Bethe ansatz equations 

We are interested in the solutions of the eigenvalue equation 

HN\'i'N)=EN\'i'N), (9) 

where is given in (8) and 

\^n)= 52 <I>(xi,X2,...,X/v)|xi,X2,...,X/v) (10) 

is the eigenfunction written in the basis that specifies the positions of the particles in the 
system, with ^{xi,X2, ■ ■ ■ ,X]\i) the coefficient for the configuration \xi ,X2,... ,xj^). Notice that 
since there is no exclusion, particle positions can coincide. 

If the positions of the particles obey xy+i > xy, I ^ j ^ N, the eigenvalue equation (9) is 
satisfied by the ansatz (10) with coefficients 

r N 

<I>(xi,X2,...,XA,) =^Ap(i)p(2)...p(/v)exp iXl^P(i)-^7 (11) 
p L j=i J 

and eigenvalue 

N 

EN{ki,k2,...,kN) = J^e{kj), (12) 

where E{k) = 1 — pe^^'^ — qe^^ is the "single particle energy," the first summation in (11) is over 
all the A'^! permutations P of the indices {\,2,. . . ,N) used to label the positions of the particles, 
and the "wave numbers" k\,k2,. ■ ■ ,kN are chosen so that \^n) satisfies (9). We see that when 
the particles are far apart, in the case being just not on the same site, they behave as if they were 
free, and the total "energy" of the system is the sum of the "energies" of single particles. When 
a pair of particles sit on the same site, Xj+\ = xj, we obtain from (9)-(12) that the amplitudes 
^p(i)P(2) - p(A') should satisfy 

Ap(l)p{2)-P(j+\)P{J)-P(N) 

where the "two-particle scattering phase" @ji is defined by 

^i0,f ^ p + qe^^^'+^^^-e'^' 

The boundary condition <I>(x2,X3, . . . ,xi +L) = <I>(xi,X2, . . . ,x/v) furnishes the additional rela- 
tion 

^P(l)P(l)-P(N) = ^'*'''*''^^P(2)P(3)--P(1)- (15) 
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Iterating relation (13) times, (15) gives us the Bethe ansatz equations for the "wave numbers" 
kj in the A'^-particle sector, 

(16) 

The solutions ,^2j • • • of these equations give through (12) the eigenvalues of (8). Notice 
that since is in general nonhermitian, the kj are in general complex numbers. 

The eigenf unctions (10) with the coefficients (11) should also be eigenf unctions of the 
translation operator T that shifts the positions of the particles to the left by one site, since 
[H, T] = 0. The eigenvalues e'^ of T are given by 



N 

E 

where we have defined the total momentum-like P by 



T\^>N)=e'''\^>M)=(jle^'^y^N), (17) 



P = 5] /cy (mod 271) = , O^^^L-1. (18) 

7=1 ^ 

With P defined above, equations (16) can be rewritten as 

^k,(L^N) ^ (_i)A^-i,ii:t... j-^,ie, ^ (-l)^-^.'^n - l^i^^- (19) 



=1 



£=1 p + qe^^''i+^f^ -e'*^' 



The learned reader will recognize in (19) the Bethe ansatz equations for the asymmetric simple 
exclusion process of particles in a lattice of L+N sites with twisted boundary conditions, the 
angle of twist being given by the total momentum P of the system. For stochastic processes, 
the relevant momentum sector is the P = sector, since the coefficients ^{x\,X2, ■ ■ ■ ,xn) have 
to be all real and positive. In this sector, the correspondence between the zero range process 
and the asymmetric simple exclusion process is exact. 



4.2 The dynamical critical exponent 

Numerical simulations together with theoretical arguments and explicit calculations indicate 
that the critical behavior of the interface is independent of the particular values of p and q 
as long as p q [13, 14, 15, 17]. It has then become usual to investigate the Bethe ansatz 
equations (16) with p = I, q = 0, since this facilitates the analysis considerably. The more 
general cases < q < p < I were investigated in [15, 36]. The p = I, q = case corresponds 
to a 2D Ising interface evolving in a finite field B > but at zero temperature. However, at 
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least for very low (but nonzero) temperatures, one has the same kind of critical behavior as 
observed at zero temperature [35]. In the totally asymmetric simple exclusion process, another 
simplification of the Bethe ansatz equations comes with the choice of the half-filled sector 
2N = L. The analogous choice for the zero range process is to consider the sector with N = L, 
which corresponds in the interface scenario to an interface with average slope 71/ 4. 

The dynamical critical exponent z, that measures the degree of anisotropy between the 
spatial and temporal correlation lengths, can be determined from the asymptotic behavior of 
the gap Ej^'' (L) of Hm through RejE"^^^ (L)} ~ L For the asymmetric exclusion process with 
arbitrary p ^ q and p = N/L, the large L asymptotic value of Ej;J\L) is given by 

Ei^\L) = -2C\p-q\ Vp(1-p)L-3/2 ±2ni\{p - q){\ -2p)\L-\ (20) 

with an exact (numerically evaluated) C = 6.509 189 . . . [36]. The dynamical critical exponent 
of the asymmetric exclusion process is then z = 3/2, indicating that it belongs to the Kardar- 
Parisi-Zhang universality class of critical behavior [18]. 

For the zero range process, the available calculations of Ej^\L) are based on the analysis 
of slightly generalized models, with non-uniform hopping rates or in which particles can hop 
together [26, 27]. The totally asymmetric zero-range process with uniform rates in the N = L 
sector was solved in [29, 30], with the result that the first gap behaves Uke 

4^)(L)~aoL-^/2-i|L-i, (21) 

with an exact (numerically evaluated) ao = 2.301345 . . ., first obtained in [14]; see also [15].* 
We thus see that, in either case, the gap behaves asymptotically as E^^ (L) ~ L^^/^, the dynam- 
ical critical exponent z = 3/2, and both processes — the asymmetric simple exclusion process 
and the asymmetric zero-range process — ^belong to the Kardar-Parisi-Zhang universality class 
of critical behavior, and so does the driven interface dynamics of the 2D stochastic Ising model 
in the regime of low temperatures and fields. 

5 Generalized particle-height model and interface dynamics 

The particle-height model mentioned in Sec. 3 can be generalized to processes in which the 
particle positions observe the constraint Xi^i —X(^s, generating a process in which particles 
move only if the next particle is far apart by at least s sites. Identifying a particle with the 
up spin state and a hole with the down spin state in the basis, the time evolution of this 

•'■Notice that the constants C in (20) and uq in (21) are related by C = 2\/2aQ. 
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restricted exclusion processes is governed by the infinitesimal generator [29, 30] 



=1 



A, (22) 



where p and ^ = 1 — p are the rates at which particles hop respectively to the right and to the 
left, 

= tl\U^- ^l) + i(l + n - ^^Vy)! (23) 

is the operator that projects out configurations in which particles are closer than by s sites, and 
= 2 (a^ ± ia^') and are the usual Pauli spin-^ matrices. 

The infinitesimal generator Hs is exactly solvable by the Bethe ansatz for any choice of the 
integer s [29, 30]. The s = \ case recovers the generator of the asymmetric simple exclusion 
process, while the generic > 1 case describes, e.g., the driven diffusion of s-mers on the lattice 
[39]. In fact, it has been shown that an arbitrary mixture of i'-mers of different sizes, diffusing 
with the same rates but ruled by evolution operators Hg with different s, can be integrated 
exactly, with the eigenspectrum depending only on the average size s of the i'-mers [29, 30]. 

For the generalized particle-height model, either with a single type of particle or with an 
admixture of particles of different sizes, the simplifying filling fraction analogous to the condi- 
tion A'^ = L in the zero-range process is given by ( 1 + s)N = L, where s is the average size of the 
particles in the system, not necessarily a semipositive integer. Clearly, the larger the average 
particle size s, the smaller the average intef ace slope N /L = l/{\+s).la this case, the spectral 
gap of the process has been found to scale like [29, 30] 

E^\L)^aoL-y' + i{^j^y-L-\ (24) 

with the same qq = 2.301345 ... as before; compare with (21). We can then predict that the 
interface dynamics obtained from the generalized particle-height model with x^+i —Xf^^s also 
belongs to the Kardar-Parisi-Zhang universality class of critical behavior. 

A most interesting thing would be to vizualize how the Ising interface configurations evolve 
in the case of a particle-height model that includes besides particles of positive sizes also par- 
ticles with negative sizes, since in this case handles and loops could develop. More generally, 
it would be of interest to find physical applications of the operator (22) for negative values of s 
or s, since they are all exactly solvable and display the same type of critical behavior. 
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6 Summary and conclusions 



We showed that it is possible to map the interface dynamics of the 2D stochastic Ising model in 
the regime of low temperatures and fields into an exactly solvable interacting particle system 
of hopping particles without exclusion. The infinitesimal generator of the process is exactly 
solvable by the Bethe ansatz, with an spectral gap in the asymmetric case scaling like L^^/^ 
with the system size. The 2D Ising interface in the presence of a driving field then grows 
and fluctuates according to the Kardar-Parisi-Zhang universality class of critical behavior. We 
remark that most studies (e.g., in the realm of nucleation dynamics) of 2D Ising interfaces are 
carried out in zero temperature, and several results were obtained in the absence of external 
fields. In our study, following [5, 35], we allow for finite temperatures and fields, as long as the 
conditions stated in section 2 are met. When p = oo (T = 0+) or B = 0, the process becomes 
symmetric and the Bethe ansatz analysis reduces to a simple spin-wave analysis. In this case 
the gap becomes ^^''(L) = 2sin^(7r/L) (independent of N as long as N < L), with asymptotic 
behavior ^^^^^(L 1) = ItI'^L^^, and the interface only fluctuates, without moving or growing, 
according to the z = 2 Edwards-Wilkinson universality class of critical behavior [37]. 

It would be desirable to explore the mapping of the Ising interface to the zero range process 
to investigate step-step correlation functions by tagged-particle methods within the context of 
exact Bethe solutions [6, 7, 9, 19, 20, 23, 38], as well as the dynamics of special configurations 
like semi-infinite strips of the minority phase ("Ising fingers") [40] using some of the ideas 
exposed here. We hope to return to these subjects soon. 
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